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SUMMARY 

The development of theoretical equations describing the 
behavior of a lumped-constant pressure-sensing system under 
transient operation is presented with experimental data that 
show agreement with the equations. A pressure -sensing system 
consisting of a tube terminating in a reservoir is investi- 
gated for the transient relation between a pressure disturbance 
at the open end of the tube and the pressure response in the 
reservoir. Design parameters are presented that can be adjusted 
to achieve a desired performance from such a system when the 
system is considered as a transfer member of a control loop. 


INTRODUCTION 

In many control applications, the controlled variables 
such as fluid flow, thrust, airspeed, torque, and temperature 
(references ! and 2) may be measured by pneumatic means. Such 
a measuring system may incorporate a tube that terminates in 
a reservoir having a pressure-sensitive element. An investi- 
gation concerned' with the transient relation between the pres- 
sure disturbance at the mouth of such a tube and the pressure 
response in the reservoir as affected by changes in system 
dimensions and gas conditions was conducted at the MCA Lewis 
laboratory and is reported herein. 

The transient behavior of linear systems may be described 
in various ways. The methods used in the present report employ: 
the differential equation describing the system, the transfer 
function (reference 3) (the Laplace transform of the response 
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divided by the Laplace transform of the disturbance), the indie ial 
response (response to a unit-step disturbance), and the frequency 
response (steady -state response to sinusoidal inputs). In synthe- 
sizing an over-all control system, the transfer function or 
frequency response has greater utility; whereas, in the study of 
a single transfer member, the transient behavior is easily under- 
stood by the use of the indicial-response method. 

Considerable research has been conducted on this problem 
in fields other than controls. References 4 and 5 report 
investigations of this problem as applied to airspeed indicators 
and altimeters. In these applications, the effect of the mass 
inertia of the flowing medium in the tube is usually so small 
as to be negligible; whereas, for controls a design that would 
permit the exclusion of the mass -inertia effect reduces the 
order of the equations and limits the flexibility of design. On 
this basis; the analyses of references 4 and 5 are not considered 
applicable to the study of pres sure -sensing systems used in 
controls . 

In the field of acoustics, sound transmission in pneumatic 
systems has been thoroughly investigated and is presented in 
references 6 and 7. Parts of these analyses are used in the 
investigation presented herein. 

Experimental data in the form of indicial responses are 
presented to substantiate the analysis. (The symbols used are 
defined in appendix A) . 


ANALYSIS 

If consideration is limited to a gas -pressure-sensing 
system (fig. l) in which the volume of the reservoir is large 
compared with the volume of the tube, the capacitance (appendix B) 
may be considered a function only of the reservoir volume and 
the gas conditions. Also, by restriction of the analysis to the 
case in which the tube length is short, the dead time (length 
of tube divided by the speed of sound) due to transportation lag 
may be neglected. As a result of these two restrictions, the 
volume flow at any time is considered uniform along the tube 
length. The "inertance" and the resistance (appendix B) are 
therefore functions only of the tube dimensions a nd the gas 
conditions. By these restrictions, the analysis is reduced to 
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an investigation of a lumped -constant system (reference 8), a ' 
system that disregards the space variables and assumes disturb- 
ances to occur instantaneous ly throughout the system with varying 
magnitudes . 


General Equations 

Differential equation describing pressure transient . - If an 
adiabatic process and small pressure changes are assumed, the 
relation between the pressure in the reservoir and the pressure 
at the mouth of the tube is described by a second-order differ- 
ential equation with essentially constant coefficients: 


JCp c + BCp c + P c = P d 


( 1 ) 


or 




p c = p d 


( 2 ) 


The definition and the derivation of the coefficients appear in 
appendix B . 

Transfer function of pres sure -sensing system. - By use of 
the Laplace transform, the transfer function for the differential 
equation is 


F(s) 


P c (s) 



1 


s + 1 


(3) 


Indiclal response . - The indicial responses in dimensionless 
form for various damping ratios (reference 9) are: 
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Case’ I: Underdamped, 0 < £ < 1 



(yrr? „ o t 



Case II: Critically damped, £ = 1 


( 4 ) 



(0 te 
0 


d) 0 t 


( 5 ) 


Case III: Overdamped, £ > 1 



These indicial responses are plotted for various damping ratios in 
figure 2 . 


Frequency response, - The equation for the frequency response 
( reference 3) is 


F(i<o) 



1 


+ 2i £ ^ 
(0 0 


+ 1 



( 7 ) 


The gain and the pihase shift as functions of frequency ratio for 
various damping ratios are shovn in figure 3. 


Variation of Coefficients 

In the design of a pressure-sensing system, consideration 
must he given to the variation of the coefficients of the differ- 
ential equation describing the system. The coefficients of equa- 
tion (2) vary in tvo distinct ways: (l) with operating pressure 
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and temperature (initial pressure and temperature of the transient), 
and (2) during a pressure transient. Because the variation is dif- 
ferent in each case, the cases must he considered separately. 

Variation of coefficients vith operating level. - The undamped 
natural frequency is a function only of the operating temperature of 
the working fluid for a given system: 

“0 "a/eT 7gR g T o ( 8 ) 

The percentage change in undamped natural frequency with percentage 
change in operating temperature is shown in figure 4. 

The damping ratio varies with viscosity, temperature, and pres- 
sure of the working fluid; 


t = -ifcL 
s 3 



( 9 ) 


Figure 5 shows the percentage change in damping ratio with percentage 
change in operating pressure. 

Because the viscosity is a function of the gas and the tempera- 
ture, a general relation showing the variation of damping ratio with 
temperature cannot be determined. The variation of damping ratio, 
however, can be shown for given gases and system dimensions. The 
variation of damping ratio with temperature for air in a system 
having dimensions such that the damping ratio is unity at standard 
temperature is shown in figure 6. 

Variation of coefficients during pressure transient, - During a 
pressure transient, for adiabatic conditions, the undamped natural 
frequency varies with the pressure as follows: 


(Or 


(0 


0, initial 


111 



( 10 ) 


The percentage change in u n d amp ed natural frequency with a percentage 
change in pressure is shown in figure 7. 
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If a linear variation of viscosity with temperature 

\i = + Kg T 

is assumed, the change in damping ratio is a function of the change 
in pressure and the initial temperature, as shown by the following 
equation: 




initial 


El + £«r 


Ki 


|T„ 



% 


JL 

W 


Zzl 

27 


(11) 


The variation of damping ratio with absolute pressure for air at 
standard initial temperature using a linear variation of viscosity 
with temperature is shown in figure 8. 


APPARATUS AND PROCEDURE 

A schematic representation of the experimental apparatus used 
to obtain transient pressure data is shown in figure 9. In order 
to obtain an approximate step pressure disturbance at the mouth 
of the tube, an air stream was interrupted by a revolving slotted 
dish. The diameter of the disk and the speed of the motor are 
such that the time of travel of a point on the disk across the 
mouth of the tube is of the order of 1/1000 second. 

The pressure response in the reservoir at the end of the tube 
was sensed by a commercial device for converting pressure to an 
electric signal. This unit, consisting of a' bellows connected to 
a strain-gage bridge, has a natural frequency of 1100 cycles per 
second. The volume flow through the tube due to the movement of 
the bellows' is negligible as compared with the volume flow due 
to the compressibility of the gas in the reservoir. 

The bridge was energized with an audio-signal generator supply- 
ing an 8000-cycle alternating-current potential. The bridge unbalance 
was fed directly into a cathode-ray oscilloscope equipped with a 
recording camera . 

The unbalance of the bridge appeared as an amplitude modulation 
of the carrier wave. This modulation was found to vary linearly with 
the pressure in the reservoir. The period of the carrier wave pro- 
vided a time base. 
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RESULTS AND DISCUSSION 

In order to determine the validity of the analysis, experi- 
mental curves of indie ial responses were obtained and compared with 
the theoretical responses calculated from equations (4) and (6). 
These curves are shown in figures 10 to 12. The damping ratios and 
the undamped natural frequencies were varied only by changing the 
tube radii. 

Close agreement between the theoretical and experimental 
responses is shown in figures 10 and 11. For these figures, the 
volume ratios (tube volume divided by reservoir volume) were 0.08 
and 0.12, respectively. 

The volume ratio was made approximately 0.50 for the data in 
figure 12. An examination indicates that the undamped natural 
frequency is lower and the damping ratio is higher than calculated. 
The volume ratio was made large in order to indicate the magnitude 
of the error involved. 

If a pneumatic system is to be used as a transfer member of a 
control system, a specific transfer function (equation (3)), indi- 
cial response (equations (4) to (6)), or frequency response 
(equation (7)) will be required. In order to attain the required 
response, definite values of undamped natural frequency and damp- 
ing ratio are specified. The dependence of the undamped natural 
frequency and the damping ratio on the system parameters are indi- 
cated by equations (8) and (9), respectively. From these equations, 
the damping ratio and the undamped natural frequency can be seen to 
vary with the operating pressure and the temperature; a system 
designed for one operating level therefore may not function properly 
at another. Because of these variations, the design must be such 
that the deviation of undamped natural frequency and damping ratio 
are within acceptable limits over the range of operating pressures 
and temperatures for which the system is designed. 

For a gas at given conditions, undamped natural frequency and 
damping ratio are f unctions of the length of the tube, the radius 
of the tube, and the reservoir volume. If practical considera- 
tions fix one of these system dimensions, the other two are uniquely 
fixed for given values of frequency and damping ratio. When two of 
the system dimensions are fixed, the third dimension is determined 
by a choice of either frequency or damping ratio. 

When design conditions are such that tube length and reservoir 
volume are at practical minimums, the undamped natural frequency 
can be increased further only by an increase in tube radius. This 
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value of tube radius may result in too low a value of damping 
ratio. In order to increase damping ratio to the desired value, 
additional resistance in the form of a restriction may be intro- 
duced into the tube. This effect is demonstrated in the example 
shown in figure 13. In this example, it is desired to have a 
damping ratio of 0.727 and an undamped natural frequency of 
571 radians per second. Practical considerations limited the 
tube length to 16.75 inches and the reservoir volume to 0.202 cubic 
inch. For the desired damping ratio, the radius must be equal to 
0.0215 inch. With this radius, the undamped natural frequency is 
279 radians per second (fig. 13, curve a). This value is lower 
than the desired frequency. Inasmuch as undamped natural fre- 
quency is a function of tube length, tube radius, and reservoir 
volume and whereas tube length and reservoir volume are at a min- 
imum, the undamped natural frequency can be increased to the 
desired value by an increase in tube radius only. This increase, 
however, decreases the damping ratio to 0.085 for a tube radius 
of ’0.044 inch (fig. 13, curve b). 

, In order to restore the damping ratio to 0.727 and to main- 
tain the un d amp ed natural frequency at 571 radians per second, 
added resistance was introduced into the tube while maintaining 
the same system dimensions (fig. 13, curve c). A method for 
increasing resistance without changing the system dimensions is 
the insertion of an orifice, a wire mesh, or other restriction 
in the tube. Inasmuch as the system dimensions are unchanged, 
the resistance is increased independently of the inertance q nH 
the capacitance. 

The damping ratio for figure 13, curve c, could not be cal- 
culated from the theoretical equations. It was approximated by 
comparing its response to curves shown in figure 2. In this way, 
the desired undamped natural frequency and the damping ratio were 
obtained, although the tube length and the reservoir volume were 
fixed. 


SUMMARY OF RESULTS 

A gas-pressure-sensing system consisting of a tube terminat- 
ing in a reservoir was considered. An analysis, restricted to a 
system having dimensions such that it can be treated as a lumped- 
constant system, was presented with experimental data. The anal- 
ysis assumed that the coefficients of the differential equation 
were constant during a pressure transient. The variation of 
these coefficients during a pressure transient can be calculated, 
however, and were presented for a particular system. 
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Agreement vas shown between the theoretical analysis and 
experimental data over the limited range of conditions investi- 
gated. Because of the large number of variables, determination 
of the limits of application of the differential equation was 
impractical. It seems unlikely that large errors would be intro- 
duced, however, if the basic restrictions were adhered to; the 
volume of the reservoir must be large compared with the volume of 
the tube, the tube length must be short, and the pressure differ- 
ences existing in the system must be small. 

Design parameters were presented, which can be adjusted to 
achieve the desired transient performance from the system when it 
is considered as a transfer member of a control. The variations 
of the coefficients of the differential equation resulting from a 
change in the operating level of the pressure or temperature were 
analyzed and presented. 

The method whereby the damping ratio of the system may be 
increased without affecting the undamped natural frequency was 
experimentally investigated. The method consisted in increasing 
the resistance to flow without changing the system dimensions. 

This method results in greater flexibility in design, inasmuch as 
a system may be designed to produce the required value of undamped 
natural frequency and any damping ratio less than that required. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, May 25, 1949. 
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APPENDIX A 


SYMBOLS 

The following symbols are used In this report: 

A area, sq ft 

a acceleration, ft/sec^ 

C capacitance, ft^/lb 

c sound velocity, ft/sec 

F(s) transfer function 

F(ito) frequency-response function 

G gain, or amplitude amplification 

g gravitational constant, ft/sec^ 

i V^l 

J inertance, lb-sec^/ft'-* 

constants 

L length of tube, ft 

mass, slugs 

Laplace transform of pressure change in reservoir 

Laplace transform of pressure disturbance 

pressure change in reservoir, Pg-p^, lb/sq ft 

pressure disturbance, P^-Pq, lb/sq ft 

pressure drop necessary to accelerate fluid in tube, 
lb/sq ft ■ ?'■■■■■ 


-0 


m 




pressure drop in tube due to flow resistance, lb/sq ft 
initial steady-state absolute pressure, lb/sq ft 
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Pi 

P2 

Q 

E 

R g 

r 

s 



U 

P 

p 0 

(0 

w o 

Subscripts : 


absolute pressure at mouth of tube, lb/sq ft 

absolute pressure in reservoir, lb/sq ft 

volume displacement in tube, cu ft 

resistance, lb-sec/ft 5 

gas constant, ft-lb/(lb) (°F) 

radius of tube, ft 

complex variable 

absolute temperature, °E 

initial steady-state absolute temperature, °B 
absolute temperature in reservoir, °R 
time, sec 

volume of reservoir, cu f t 
ratio of specific heats 
damping constant, sec" 1 
damping ratio 
phase shift, radians 
absolute viscosity, lb-sec/sq ft 
density in tube, slugs/cu ft 

initial steady-state density in tube, slugs/cu ft 

input frequency, radians/sec 

undamped natural frequency, radians/sec 


a,b operating levels 

Dots above the symbols represent derivatives with respect to time. 
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APPENDIX B 


DERIVATION OF THEORETICAL EQUATIONS 

The pressure-sensing system shown in figure 1 can he treated 
as a lumped-constant system and described by an ordinary differ- 
ential equation using the following restrictions and method of 
analysis. 


General Equation 


Restrictions. - 


(1) The tube length must be sufficiently short so that the dead 

time — can be neglected, 
c 

(2) The dimensions of the reservoir must be such that the 
pressure throughout the . reservoir may be considered uniform at 
any time. 

(3) The volume of the tube must be small compared with the 
volume of the reservoir. On this basis, the assumption can be 
made that the volume flow is uniform throughout the length of the 
tube and is the result of the compressibility of the gas in the 
reservoir. 

Definition of capacitance C (references 1, 6, 7, and 10). - 



The volume flow into the reservoir is equal to the decrease 
in volume of the gas originally in the reservoir. 


therefore 


dQ = - dV 



When adiabatic compression is assumed in the reservoir 
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^ - H 

Y^dp 2 + p 2 rV 7-1 dV - = 0 


dV _ -Vy _ -V 

^2 " p 2 rv r - 1 = V 


Therefore 



The value of C varies inversely with the absolute pressure in the 
reservoir. For small pressure changes, C is assumed to be con- 
stant. 


Definition of resistance R. - When the resistance is assumed 
to exist only within the tube, the resistance to flow is defined as 
(references 1, 6, 7, and 10) 



From the Hagen-Poiseuille law 


Therefore 


Q = 


«P r r 

8nL 


8ixL 

irr 


The value of E varies with the absolute viscosity, which is 
assumed essentially constant for small changes of pressure tem- 
perature. The Hagen-Poiseuille law is valid only for laminar flow. 
For small Reynolds numbers, R is constant. 

Definition of lnertance J. - If the motion of the fluid in 
the system is assumed to occur only in the tube, the inertance is 
defined as (references 6, 7, and 10) 
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m = pLA 


The acceleration of the fluid in the tube is 




Then 


Therefore 




j = £k = jeL 

a ^2 

itr 

The inertance is proportional to the density. For small pres- 
sure and temperature changes, the inertance is essentially constant. 

Differential equation. - The following equation can be obtained 
from a pressure balance across the tube: 

P J + p r + p 3 * P 1 


Then 


p j + P r + P 2 - P 0 ?! - P 0 


p j + p r + p c - Pa 


P J 

= JQ 

p r 

= BQ 

dp 2 

C 

dp 2 

■If- 

P 2 

= K 4 


By definition 
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when 



Q = 0 


Therefore 

p 2 “ p 0 = C = P c 

and 

JQ + EQ + | = p d 


JC P c + EC p c + p c = p d (Bl) 

The general nature of this equation is so well known tha t a 
detailed discussion is unnecessary. By definition 


(0 n = — — (undamped natural frequency) 
u VJC 


5 = 2 j (damping constant) 

£ = ~ = | (damping ratio) 


Equation (Bl) reduces to 


P + 2 4- p + 


o P. 
00 * c 
0 


0 ) 0 p c 


p c = p d 


(B2) 


Transfer function. - , By use of the Laplace transformation, the 
transfer function is (references 3 and 8) 


F(s) = 


F c( s ) 1 

p d( s ) s 2 „ £ 

+ 2 rf- S + 1 


(B3) 


0) 2 
0 


Cl) 
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Indie ial responses. - For a pressure disturbance p d equal to 

a step input, the solutions of equation (B2) in dimensionless form 
are (reference 9): 

Case 1: Underdamped, 0 5 £ < Is 


p c 

— = 1 - * - ■ 

Pci * A _ >2 


-i £ 


vrr? cos ? 


(B4) 


Case II: Critically damped, £ = 1: 


— = 1 - e 

Pd 

Case III: Overdamped, £ > 1: 


-<*>nt 


- 0)r>te 




(B5) 


h . 1 - e' t%t 

P. 


sinh A^£ 2 - 1 w Q t + cosh ^£ 2 - 1 a) Q t 


a/£ 


(B6)" 


Frequency response. - The frequency response is obtained from 
the transfer function (equation (B3 ) ) by substituting 



Variation of Coefficients of 
Differential Equation 
Variation with operating level. - 


1. Undamped natural frequency 
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When the working fluid is assumed to he air, and the viscosity 
is assumed to vary linearly with temperature and is independent of 
pressure (reference 11) 

\i m 88.9 X 10" 9 + 0.548 X lO -9 T Q 

4(88.9 + 0.548 T Q ) 10' 9 
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For a change in operating pressure. 


£a 

r 


'0,b 


0,a 


If the system dimensions are such that £ = 1 at standard tempera- 
ture, the magnitude with operating temperature is 


£ = (1.176 X 10" 4 ) (88.9 + 0.548 T Q ) ^ 


Variation during pressure transient. - In order to simplify the 
analysis, the assumption is made that only small pressure differences 
exist in the system at any time. 

1. Undamped natural frequency 


_ J2_ 

' ^g T 2 


U ) 0 




If an adiabatic process is assumed during the transient 


Zli 

®o Jh ) 27 

W 0, initial \ p 0/ 


(BIO) 


2. The damping ratio is 
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assuming 


n - Ki + Kg T 2 
4(K 1+K2 T 2 ) 


t- 


P 2 r ^ 


*7 ® R 6 T 2 


during a transient 


t 


K 1 + *2 T 2 / p 0 \ /^ T 2\ 2 
tJ 


£ initial *1 + *2 T 0 \ P 2, 

assuming adiabatic conditions during a transient 


i 


initial 



For air. 


_J 

►initial 



-(2± l) 

'P 2 s 


(Bll) 
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0 1 2 3 4 5 6 7 8 


Dimensionless time, o> 0 t 

Figure 2. - Dimensionless indicial response of second-order linear 
differential equation with constant coefficients. 
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400 440 480 520 560 600 640 


Operating temperature, Tq, °R 

Figure 6. - Variation of damping ratio of pressure -sensing system 
with operating temperature of various constant values of operating 
pressure* Operating fluid, air; damping ratio, 1 at standard 
temperature* 



Change In undamped natural 

Change in damping ratio, £ , percent frequency, co Q , percent 
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-16 -12 -3 -4 0 4 8 12 16 


Change In absolute pressure, p 2 , percent 

Figure 7. - Variation of undamped natural frequency of pressure -sensing system with 
change in absolute pressure during pressure transient. 



Change in absolute pressure, p 2 , percent 


Figure 8. - Variation of damping ratio of pressure -sensing system with change of 
absolute pressure during pressure transient. Operating fluid, air; initial 
steady-state absolute temperature, 59° F. 
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Figure 9. - Experimental apparatus used to obtain transient pressure 
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p d/°d ‘asuodeau 
IBfOIPlII 


p d/°d ‘esuodeaa 
I^TOTPUI 


Figure 10. - Indicial response of system with calculated damping ratio of 1.35. Undamped 
natural frequency, 227 radians per second; initial pressure, 29.32 inches mercury absolute 
pressure disturbance, 1.00 inch mercury; initial temperature, 525° R; tube radius, 0.0175 
inch; tube length, 16.75 inches; reservoir volume, 0.202 cubic inch; volume ratio, 0.08. 
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Pd/°d 'asuodsaa 

IBfOTpui 


p d/°d ‘asuodsaa 

I^TOTPUi 


Figure 11, - Indicial response of system with calculated damping ratio of 0*727. Undamped 
natural frequency, 279 radians per second; initial pressure, 29.32 inches mercury absolute 
pressure disturbance, 1.00 inch mercury; initial temperature, 525° R; tube radius, 0.0215 
inch; tube length, 16.75 inches; reservoir volume, 0.202 cubic inch; volume ratio, 0.12. 
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Pd/^d * dsuods 9*i 
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^d/°d ‘asuodsaa tejofpui 


(b) Plotted values. 

Figure 12. - Indicial response of system with calculated damping ratio of 0.085. Undamped 
natural frequency, 571 radians per second; initial pressure, 29.32 inches mercury absolute 
pressure disturbance, 1.00 inch mercury; initial temperature, 525° R; tube radius, 0.044 
inch; tube length, 16.75 inches; reservoir volume, 0.202 cubic inch; volume ratio, 0.50. 
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Pd/°d ‘esuodsea x B T°TP u I 


NACA-Langley - 12-20-49 - 1000 


gure 13. - Theoretical and experimental indicial responses illustrating effect of 
inserting additional resistance in tube. Initial pressure, 29.32 inches mercury 
absolute; pressure disturbance, 1.00 inch mercury; initial temperature, 525° R; 
tube length, 16*75 inches; reservoir volume, 0.202 cubic inch. 




